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I. INTRODUCTION 



In the last few years a lot of work has been investigated to study the behavior of vector 
mesons in a medium with finite baryonic density. The basic motivation was to find a sign of 
chiral symmetry restauration in heavy-ion collision experiments, when studying the dilepton 
spectra which correspond to the vector mesons. Indeed, the CERES experiments for S- 
Au and Pb-Au collisions show a novel feature when compared to proton-nucleus collisions, 
namely an enhancement of the dilepton yield for invariant masses somewhat below the 
vacuum mass of the p meson |]-f§. Some years ago it was argued by Brown and Rho that 
such an enhancement might be due to the restoration of chiral symmetry || . In their model 
they assumed that the masses of the vector mesons should scale with the quark condensate, 
i.e. drop with rising baryonic density. If this is true, the p peak in the dilepton spectrum 
would be shifted to lower invariant mass. This might be an explanation for the observed 
enhancement of the dilepton yield in that region However, the idea of chiral symmetry 

restauration alone without additional model assumptions does not provide a unique picture. 
Other scenarios predict a rising p mass based on the effect that the p becomes degenerate 
with its chiral partner, the al meson ||. 

Besides the problem what the consequences of chiral symmetry restauration might be 
there is still the possibility that the experimental finding of the enhancement might also 
be explained by conventional hadronic degrees of freedom. To clarify that issue hadronic 
models for the in-medium behavior of vector mesons were developed by various groups, see 
e.g. p|-|18[. Some of them predict a large peak broadening of the p meson or even distinct 
new peak structures. It was found that the enhancement in the dilepton yield might also be 
explained within a purely hadronic scenario, if a lot of strength is shifted to lower invariant 
mass [|l4tp!9| . So far one is not in a position to confirm or rule out such hadronic models 
by experimental data. Therefore, it is of interest to find additional model independent 
consistency checks which should be obeyed by arbitrary hadronic models describing vector 
mesons and their in-medium behavior. Such a consistency check is provided by the QCD 
sum rule approach. 

Originally, QCD sum rules were developed for vacuum processes not as a consistency 
check for hadronic models, but as an alternative to them, i.e. to deduce model independent 
information about hadrons from the underlying quark and gluon degrees of freedom (see 
e.g. [pOH 22[| ) . An important ingredient for the description of vector mesons by QCD sum 



rules is the assumption that the spectral functions of these resonances can be reliably ap- 
proximated by ^-functions (narrow width approximation). In this way, the experimentally 
found vector meson masses can be reproduced reasonably well. Of course, it is supported 
by experiments that e.g. the respective width of the spectral function of p and u meson is 
small as compared to the mass of the meson. However, it is important to realize that the 
narrow width approximation is not a result, but an ingredient of the traditional QCD sum 
rule approach. 

In the last few years QCD sum rules were also developed for in-medium situations, 
i.e. for hadronic matter at finite temperature (see e.g. |23| , |2~4}| ) or finite baryonic density 
(e.g. |25|-^9|). For p and u mesons it was found that their masses decrease with increasing 



temperature and/or density, if the narrow width approximation is used for the parametriza- 
tion of the respective spectral function. In contrast to the vacuum case, there is, however, no 
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experimental support that the narrow width approximation is a reasonable assumption for 
the in-medium case. Unfortunately, this approximation crucially influences the QCD sum 
rule prediction for a possible mass shift in a nuclear surrounding. If a spectral function with 
an appropriately chosen large width is used in the QCD sum rule approach at finite density, 
one could get an unshifted meson mass, in contrast to the finding utilizing the narrow width 



approximation. This was first discussed in |15[ using a specific hadronic model and later 
systematically studied in |3"0fl . 

This shows that QCD sum rules provide no model independent prediction about a pos- 
sible mass shift of vector mesons in nuclear medium. Only together with some additional 
assumptions (e.g. about the width of the respective vector meson) a statement about the 
density dependence of the masses of the vector mesons can be deduced from the sum rule 
analysis. Nevertheless, once a hadronic model for vector mesons has been chosen the sum 
rules can be used as a consistency check for this model. We believe it to be important to 
have such consistency checks, since it is not clear a priori if a hadronic model — e.g. with 
coupling constants derived from vacuum processes — yields a correct description of the 
in-medium behavior. 

In most of the studies on QCD sum rules only the vector mesons which are at rest 
with respect to the medium were considered. On the other hand, since Lorentz invariance 
is broken, the behavior of vector mesons clearly depends on their velocity with respect to 
the surrounding. Indeed, some of the hadronic models mentioned above yield very differ- 
ent spectral functions for different three-momenta of the respective vector meson and for 
different polarizations (e.g. ||17|| ). Only recently, the influence of finite three-momentum on 
the QCD sum rule prediction for vector meson masses (within the classical narrow width 
approximation) was explored [pE]. Also at finite three-momentum, of course, the narrow 
width approximation is an assumption which may be justified or not, but in any case it is 
not a model-independent statement. 

In this work we will derive QCD sum rules for p and u mesons for arbitrary three- 
momentum of the respective vector meson with respect to the nuclear medium. The purpose 
is to provide a consistency check for hadronic models. This is in contrast to the traditional 
QCD sum rule approach which aims at a prediction for the vector meson mass assuming that 
the width of the vector meson is negligibly small (cf. especially concerning the extension 
to non- vanishing three-momentum within the traditional QCD sum rule approach). In 
the application of the traditional QCD sum rule approach to nuclear matter the attention 
was focused on the utilization of the sum rule within the narrow width approximation 



25| , |26|p8| , |29[| , rather than on a detailed discussion of the derivation of the sum rule and of 



the calculation of the various condensates which contribute. In the present article we try to 
bridge this gap. 

In the next section we introduce the basic quantity of interest, namely the current- 
current correlator and present a dispersion relation which connects the calculations for this 
correlator using hadronic degrees of freedom on the one hand side and quarks and gluons on 
the other. In section |IT| we sketch the method of operator product expansion and calculate 
the current-current correlator within that framework. In section |Iy| we present a QCD sum 
rule derived from the dispersion relation mentioned above. To get more insight into the 
various contributions calculated in section |T| we discuss in section [V| an approximation 
linear in the nuclear density. In section [VI| we discuss the various approximations which 
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have led to the results presented in the preceding sections. Finally we summarize our results 
in section VII. 



II. THE CURRENT-CURRENT CORRELATOR 

The quantity we study in the following is the covariant time ordered current-current 
correlator 

IV(g) = i Jd A xe^(Tj^x)jM) ■ (2-1) 

Here is an electromagnetic current with the isospin quantum number of the respective 
vector meson, 

3n = o T dl^d) (2-2) 



where — is for the p meson and + for the uj. The current-current correlator enters e.g. the 
cross section of e + e~ — > hadrons (see e.g. fI5(l). Within a simple vector meson dominance 
(VMD) picture the current ( |2.2| ) can be identified with the vector meson which carries the 
respective isospin, e.g. for the p meson |fl]]: 

2 

. VMD m p , n q \ 

= — (2-3) 
9p 

where p M denotes the p meson field amplitude and g p the coupling of the p meson to pions. 
Therefore, within simple VMD the current-current correlator is proportional to the propaga- 
tor of the respective vector meson. Speaking more generally, i.e. without referring to simple 
VMD, the vector meson propagator is closely related to the current-current correlator. 

The expectation value in ( |2.1| ) is taken with respect to the surrounding medium. We 
study here a (isospin neutral) homogeneous equilibrated medium with finite nuclear density 
and vanishing temperature. In the medium Lorentz invariance is broken. All the formulae 
we will present in the following refer to the Lorentz frame where the medium is at rest, 
i.e. where the spatial components of the baryonic current vanish. 



The current-current correlator can be decomposed in the following way 31]: 



IV(g) = U T (q)T^(q) + IL L (q) L^(q) , (2.4) 

where we have introduced two independent projectors L^ v {q) and T^ v {q) which both satisfy 
current conservation g M L Aljy (g) = q^T^q) = and add up to 

T lv {q)+L lv {q)=g lv -^-. (2.5) 

The tensors T and L are transverse and longitudinal with respect to three-momentum q, 
respectively. T is given by 



,/i = 0orz/ = 
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while L can be deduced from fl2.5|) . The scalar functions Ht,l can be obtained by 

^(g)T^(g) = i^ + | 

and 



n T (g 2 , g 2 ) = ~iV(g)T<-(g) = ~ f + f^rO (2.7) 



II L (g 2 , g 2 ) = UMlTiq) = ~^U 00 . (2.8) 

To get the respective last equality in the last two equations, use is made of the fact that 
U lll/ (q) is a conserved quantity, i.e. transverse with respect to the /bur-momentum q. Note 
that 11^ and IT^ depend only on the invariant mass squared, q 2 , and on the three- momentum 
squared, q 2 . The latter property is due to the remaining 0(3) symmetry of the equilibrated 
system. 

At vanishing temperature the scalar functions Ut, IT^ deduced from the time ordered 
current- current correlator ( |2.1| ) can be related to the commutator (spectral function) 

AM := -~ foe^ikj^x), j„(0)]> (2.9) 
in the following way (see e.g. p2|): 

Imn TjL (g 2 , q 2 ) = -sgn(g )^ T , L (g , q 2 ) (2.10) 

where At and Al are deduced from A^ u analogously to Q2.7|) and fl2.8|) , respectively. At first 
sight, it seems that ImiT Tji does not only depend on q 2 and q 2 as claimed above, but also 
on the sign of go- However, it is easy to check from the definition ( |2.9| ) and the symmetry 
properties of the system under consideration that At and Al are antisymmetric with respect 
to the transformation go ~~ > — <Zo- Therefore we define 

A^ >L (q 2 , f) := sgn(g )^l Tii (g , g* 2 ) . (2.11) 

Inserting this relation in ( J2.10Q , 



Imn T , L (g 2 , g- 2 ) = -^, L (g 2 , g* 2 ) , (2.12) 

it becomes obvious that the dependence on the sign of go is only apparent. 

For g 2 the current- current correlator ([2.1|) can be calculated using Wilson's operator 
product expansion (OPE) |33[] for quark and gluonic degrees of freedom []2D| (for in-medium 
calculations see e.g. p3| , p5| , p6| ) . In the following we shall call the result of that calculation 
n!^ E . On the other hand, a hadronic model (e.g. for vector mesons |9]- |r5| , |r7| | using one or 
the other form of VMD) can give an expression for the current-current correlator valid in 
the time like region q 2 > 0. We denote the result of the hadronic model by IT^ £ D . A second 
representation in the space like region which has to match can be obtained from II^ D 

by utilizing a twice subtracted dispersion relation. We find 

IW, g- 2 ) = n T , L (0, g"* 2 ) + c T , L (g- 2 )g 2 - Q - ds 7'^' \ . ,, 2 (2.13) 

7r J [s — q 2 — ie){s + it 1 ) 1 
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with the subtraction constant 



2 dU T , L (k 2 ,q 2 ) 

CT,L(g ) = 



d(k 2 



(2.14) 

k 2 =0 



As we shall see below *4.p£ D (g) diverges linearly with q 2 . Therefore we have used above a 
twice subtracted dispersion relation. In the space like region for Q 2 := —q 2 >0we get the 
following connection between the current- current correlator calculated from OPE on the one 
hand side and from a hadronic model on the other: 

D 4 °r TmTT HAD f<i a 2 ) 

ngW, f) = n?r<o, f ) - crAW + ^ /* ( , + <y ? fe fcV*? (™) 

where we have used ( J2.12[ ) to express the spectral function in terms of the imaginary part 

of n. 

In the next section we shall elaborate on the calculation of the l.h.s. of ( |2 . 1 5| ) . 



III. OPERATOR PRODUCT EXPANSION 

Within the method of operator product expansion we have to calculate the current- 
current correlator ( |2.1| ) for large space like momenta Q 2 = —q 2 . Here the relevant length 
scale for the x-integration given by the inverse of \j\fQ 1 is small. This defines the hard 
scale in our problem. Suppose that the distance x is much smaller than the typical length 
of the system (soft scale). The latter might be characterized e.g. by the average particle 
distance in the medium or 1/ Aqcd as the scale where non perturbative effects appear. If x 
is that small it is reasonable to assume that a product of local operators A and B, i.e. 

A(x)B(0), (3.1) 

should look like a local operator, since the system cannot resolve such small distances x. 
Thus we find 

A(x)B(0)^C n (x)a, (3.2) 

n 

where C n denotes c-number functions (Wilson coefficients) and O n local operators. The only 
dependence on the system under consideration enters via the respective matrix elements of 
the local operators O n . Thus, the dependence on the soft scale is entirely given by the local 
operators. On the other hand, the Wilson coefficients C n can be calculated independently 
from the system under consideration. Since the operators O n are local, the dependence on 
the hard scale (here given by the x-dependence) enters only the Wilson coefficients. Thus 
we have achieved a separation of the hard from the soft scale. 

In our case, the operators A and B are the currents j M and j v , respectively. The Fourier 
transformation which appears in ( |2.1| ) does not touch the local operators O n , but only 
changes the x-dependence of C n into a g-dependence. From the line of arguments one can 
already guess that the Wilson coefficients finally yield a power series in 1/Q 2 (corrected 
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by logarithms from renormalization). The expectation values of the local operators O n 
(condensates) show up as coefficients of that series. On dimensional grounds it is obvious that 
the higher the dimension (in terms of masses) of a condensate is, the more it is suppressed 
by powers of l/Q 2 - 

In the following we will consider condensates up to dimension 6. In vacuum only scalar 
condensates contribute. For the case at hand, however, the condensates might also carry 
spin, since Lorentz invariance is broken. It is common practice to classify the condensates by 
their dimensionality d and their twist r. The latter is defined as the difference of dimension 
d and spin s, i.e. t = d — s. We decompose the current- current correlator ( J2.1|) in the 
following way: 

n OPE w nscalar + ^=4,^=2 + ^=6,^=2 + (3 3) 

where we have neglected contributions from higher dimensional condensates. We will discuss 
the various contributions separately in the following subsections. 



A. Scalar condensates 



The contribution of the scalar condensates to the current- current correlator for the system 
with finite nuclear density is formally identical to the vacuum case. The only difference is 
that in the former case the expectation value is taken with respect to the medium. The 
latter case was already discussed in the original paper by Shifman et al. [20]. The result is 
(cf. e.g. |23| for details) 



n 



scalar 



(?) 



Q 2 R scal! " r (Q 



(3.4) 



with 



i? scalar (Q 2 ) 



sir* 
na 



a 



—G 2 



2Q e 



na 



7r J \n< j 2Q 4 \ / 24Q 4 \ 7T 

(u^, l5 X a u =F d 1 ^ 5 \ a d)(uYl 5 \ a u =F dYl 5 X a d) 



(u^\ a u + dj,\ a d) £ ^A a g 

\ q=u,d,s I 



(3.5) 



where we have neglected contributions quadratic in the light current quark mass m q as well 
as differences in the masses of up and down quark. This is reasonable, since the hard scale 
is typically of the order of 1 GeV (the order of magnitude of the considered vector 
meson), while the masses of up and down quarks are of the order of a few MeV. Again the 
— sign refers to the p meson and the + to the u. To simplify (|3.5| ) we assume that the quark 
condensates of up and down quarks are approximately the same. Furthermore, we replace 
the four-quark condensates by products of two-quark condensates. Since it is not clear how 
accurate the assumption of ground state saturation (Hartree approximation) is, we multiply 
the result with a (still to be determined) factor k. We end up with 
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.scalar N _ 1 (, , M 1 J ^ 1 i m «/^\ , 1 1127ra * /^\2 



«— W „ __ ^ + ) m ^ j + _| m + _ ^ _ _^ . (3.6) 

Note that there is no difference between p and u any more, since there are no terms like 
(ud) in an isospin neutral medium. 

Of course, the crucial question is, how to evaluate the expectation values with respect 
to the nuclear medium. If the density is small, it is reasonable to approximate the medium 
by a Fermi gas of free nucleons, i.e., 

ifci<fc F 

where (O)o denotes the vacuum expectation value of an arbitrary operator O, hp the Fermi 



momentum, Ek = yfn% + k 2 the energy of a nucleon, and \N(k)) a single (isospin averaged) 
nucleon state with momentum k normalized according to 

(N(k)\N(k')) = (27r) 3 2E k 5{k — k') . (3.8) 



We will use the approximation ( ft.7| ) throughout this work and comment on it in section [VI. 

If O is a scalar operator, the expectation value (N(k)\0\N(k)) is a scalar as well and 
therefore independent of the momentum of the nucleon. Thus we get 

/d k 
fo , 3o;? ■ (3.9) 

\k\<kp 

For the evaluation of the condensates in ( |3.6| ) we need to know the expectation values of 
the quark and gluon condensate with respect to single nucleon states. The former can be 
related to the nucleon sigma term pSflfl 

(iV(0)|gg|iV(0)) = ^^, (3.10) 
m q 

while the latter can be calculated from the trace anomaly of QCD, 

(N(0) \^G 2 \ iV(0)) = -fm N mP . (3.11) 

Here denotes the nucleon sigma term and m$ the nucleon mass in the chiral limit. 

Finally, we have to make a choice for the parameter k which parametrized the deviation 
of the four-quark condensate from the product of two-quark condensates. Even for the 
vacuum case, the question about the value for k is not settled yet (cf. e.g. [p5| , p^ , p^ , p!5| ,p^1 ) . 



In addition, k might depend on the nuclear density. For simplicity, we take in the following 
the vacuum value for k also for arbitrary finite densities, but note that this introduces an 



1 Note that the normalization of the nucleon state in [25] is different from ours. 
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uncertainty into the evaluation of the OPE. All the parameters not specified so far are taken 
from []15! and listed in table |. We discuss this choice in section |V|. 

Using ( j2 . 7| , p . 8| , |3~4| . |3 . 6|]3~9^ - |3 . 1 1| ) we are able to calculate the contribution of the scalar 
condensates to the l.h.s. of |2.15|) . In the next subsection we discuss the contribution of the 
twist-2 condensates with d=4. 



B. Twist-2 spin-2 condensates 

In vacuum only scalar condensates contribute to the current- current correlator since there 
is no Lorentz vector which can account for the spin of a non-scalar condensate. Contrary to 
the vacuum case, in a nuclear medium the baryonic current can yield the spin. Using the 
approximation ( |3.7| ) we find for spin-2 condensates 

( ^ } ^ 4 I (270 3 k 2E k {N{%)P ^ N{%)) - (3 ' 12) 



|&|<ft 



F 



In this approximation the four-momentum of the nucleon accounts for the spin of the con- 
densate. Thus we get 

(N(k)\0^\N(k)) ~ hpK - 9 -fm 2 N =: S^(k) . (3.13) 

Note that the non-scalar operators are traceless with respect to the Lorentz indices. 

Expectation values of twist-2 condensates with respect to single nucleon states as they 
appear in ( |3.12| ) are thoroughly studied in deep inelastic scattering (DIS), albeit for a some- 
what different kinematical situation. We can utilize the results obtained there for our case 



at hand — as already pointed out in p5| , |26| , p9| . Therefore, we will not give a detailed cal- 



culation for the contributions of these condensates to the current-current correlator, but 
instead present a recipe how to deduce the necessary information from the DIS calculations 



of M|. 



The twist-2 operators of dimensionality 4 which contribute are given by 

ST i{u ltl D u u + d^D v d) (3.14) 

and 

STG\G KV . (3.15) 

Here ST denotes an operator producing an expression which is symmetric and traceless with 
respect to the Lorentz indices /i and v. D u is the covariant derivative and G KV the gluonic 
field strength tensor. In principle, composite operators mix under the renormalization group, 
if they have the same quantum numbers ||37|| . To study that mixing we have to decompose 
( p.!4j ) in a flavor singlet part which mixes with the gluonic operator (|3.15|) and a flavor 



non-singlet part which does not mix. For the energy region of interest, i.e. roughly about 
the masses of p and u meson, we have to deal with three active flavors. Therefore, we 
decompose ( |3.14j ) according to 
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u + d= - \2(u + d + s) + (u + d-2s)} , (3.16) 

where u is an abbreviation for u^^DyU, etc. Renormalization group mixing applies to (u + 
d + s) and G, i.e. schematically 

1 2 
w + d^-[2(w + d + s + G) + (w + d- 2s)] = u + d + -G . (3.17) 



The contribution of the twist-2 spin-2 condensates to the current-current correlator ( |2.1|) 
can be written as 

\k\<k F * 

with the twist-2 spin-2 contribution to the forward scattering amplitude between a nucleon 
and a quark current, 

Tj )s=2 ' T=2 (g, k) := i |Ae^(iV(fc)|T^(x)7^(x)^(0)7^(0)|iV(fc)) s=2>T=2 . (3.19) 



The latter quantity is calculated in [34| for the DIS case. Note that the factor 1/2 in the 
definition of the current (|2.2j ) — which enters (|2.1| ) quadratically — is not contained in the 



definition of the forward scattering amplitude, but is given in ( 3.18 ) explicitly. 

The strategy to utilize the DIS results for the forward scattering amplitude (|3.19| ) is the 
following: We make a general ansatz for this amplitude, specify it to the DIS case, match 
it with the calculations of [34], and determine in this way the unknown quantities of the 
general ansatz. 

The Lorentz structure of the forward scattering amplitude (|3.19 ) must be built up from 
the tensor g^, the four-momentum q of the quark current, and the tensor S^ u (k), defined 
in ( |3.13|) . In addition, the amplitude has to obey current conservation. Finally, it must 
be symmetric with respect to an exchange of the Lorentz indices. A general ansatz which 
fulfills all requirements is 

r W-=2,r=2 ((7j k ) = B 1 [q A S, v {k) - q 2 q,q a S ua {k) ~ q 2 q„q a S , a (k) + g , u q 2 q a q P S a p{k)] 

+ B 2 {q^q v - g^q 2 ) q a q f3 S a ^k) (3.20) 

with so far arbitrary coefficients B\ and B 2 which might depend on q 2 and q ■ k. 

The kinematical situation of DIS is such that both — q 2 and q ■ k are large with the 
Bjorken variable x = —q 2 /(2q ■ k) fixed. In this limit only the k^k v term of S fiu (k) has to 
be taken into account .0 We get 



2 Note that this is not true for the case we are actually interested in. Because the three-momentum 
q might be small, we also have to take into account the g^ v term. However, for the determination 
of the coefficients B\ and B 2 this does not matter. 
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T Ws= 2 ,.=2DIS (g; k) = _ q 2 {g . k) 2 {Bi ^ + ^ £ ^ 

with the tensors [3~3| 



and 



. + k u q^ 

(q ■ k) 2 q ■ k 

By comparison of (|3.21|) with equation (2.4) of [34] we find 

4 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



and 



Efii At 



(3.25) 



with the process independent coefficient function C l rn and the {n — l)th moment A l n of the 
distribution of the parton i in the nucleon. Expressions for the former can be found e.g. in 
[ 35| including a s corrections.^ The latter is given by [^9] 



(3.26) 



for quarks and 



A G n = 2 Jdxx n ~ l G(x,n 2 



(3.27) 



for gluons. The parton distributions xp, ip, and G at the renormalization scale fi 2 are 
parametrized in P6 |. 

Now we collect all the obtained information to get the contribution of the twist-2 spin-2 
condensates to the current-current correlator: 



nJf< T=2 (g) 



d 3 k 



\k\<k f 



(27T)3 2E fc 



T^ +d)s=2 ' T=2 (q,k) 



(3.28) 



3 Note that our notation (basically adopted from |34j]) is somewhat different from the one of [35]: 
Our coefficient function Cr, n is identical to C^j of [f35| with n = N . The longitudinal coefficient 
functions differ by a factor 2: C 3 L n = •, where again the former is our coefficient function and 



the latter the one of [35 
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with 



rp(u+d) S = 2,T = 2 
flU 



Q 

X 

4 



[q 2 S^(k) - q^q a S ua (k) - q„q a S, a (k) + g^q a q P S af5 {k)] 



r 1 '! AU+d 
L/ 2,2 /1 2 



g°2,2 /1 2 



(QnQv ~ 9^q 2 ) q' 



G 



(3.29) 



where we have taken into account that the gluonic contribution enters with a factor 2/3 
according to ( |3.17| ). The coefficient functions and the moments of the parton distributions 
are listed in table ||. We note in passing that the momentum integrations in ( |3.28|) can be 
performed analytically. Since it is not illuminating to present the lengthy result of these 
integrations we stick to the compact form given by (|3.28|) and fl3.29|). 



C. Twist-2 spin-4 condensates 

The twist-2 spin-4 condensates can be treated in the same way as the twist-2 spin-2 
condensates. We use the approximation (|3.7|) to find 



(CW>«4 / { J^ 2Ek (N(k)\0, UKX \N(k)) ■ (3.30) 

|k|<fcp 

For the expectation value with respect to a single nucleon state we get the decomposition 
(N(k)\0^ x \N(k)) 

~ kf,k u k K k x - -{k^k u g K xm 2 N + 5 permutations) + —{g^g K \m% + 2 permutations) 
=: S^k) . (3.31) 

The relevant operators are 

ST i{u lpL D u D K D x u + d ltl D u D K D x d) (3.32) 

and 

STG a ^D u D K G aX . (3.33) 
By comparison with the DIS calculations we find 

n^= 2 (g) = / T^^ T !^ +d)s ^ 2 ^ k ) ( 3 - 34 ) 

\k\<k F 

with 
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1 6 

T^ +d)s=A ' T=2 (q, k) = - — (g^gV - gffqrf - g v a q 2 q^ + V) ffV<W*) 

X IU 2| 4^ -|- g<-'2,4 /1 4 

16 

—^(q^qu - 9^q 2 ) q a q^q 1 q 5 S aPl5 (k) 



x 



(Cl,4 — Cl A )A\ +d + - (Cx (4 — C^A^j , (3.35) 

where the tensor S a p~,s(k) is defined in (|3.31|) and the coefficient functions and parton dis- 
tribution moments are listed in table [TT[ 

D. Twist-4 spin-2 condensates 

Finally we turn to the higher twist condensates. For our calculation up to dimensionality 
d=6 we need condensates with twist r=4 and spin s=2 (besides the higher twist condensates 
which are scalar and already discussed in subsection [III A|) . These condensates are [j3S| -|4l"l 

0\ v := ST g 2 ^(u W \ a u =F d W \ a d)(uw 5 \ a u =F d lul5 X a d) , (3.36) 
Ol v :=STg 2 - A (u 1 ,\ a u + d 1 ,\ a d) £ <lluX a q, (3.37) 

q=u,d,s 

0% := STig ~(u {D„ G^j^u + d{D„ G ua }j a l5 d) , (3.38) 

STg-^u [D M , G va ] + d[D„, G va \ j a d) , (3.39) 

and 

ST-fauuD^DvU + madD^Dvd) , (3.40) 

where the — sign in (|3.36|) corresponds to the p and the + sign to the u meson. For an 
unknown reason the condensate given in ( ^.39[ ) actually does not contribute to the current- 



current correlator [|39|; we have listed it here for the sake of completeness, only. In the 
following we will neglect the condensate ( |3.40| ), since it is proportional to the very small 
light quark masses and therefore suppressed. Additionally, we will neglect the contribution 



of the strange quarks to the nucleon expectation values of the operators given above JRJ . 

As in the previous subsections we use the Fermi gas approximation which for the case of 
spin-2 condensates we have already given in (|3.12[) . We can also use the Lorentz decompo- 
sition from (|3.13|) : 

(N(k)\Ol u \N(k))=:A^S, u (k) (3.41) 
where the index j = 1, 2, g denotes the twist-4 operators specified in (|3.36| - |3.38|) .n 



J The coefficients A> should not be confused with the parton distribution moments A\. 
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In j|0] it is thoroughly discussed how the expectation values of these condensates with 
respect to single nucleon states can be extracted from DIS experiments. To determine 
the contribution of these condensates to the current- current correlator we follow the same 



strategy as described in subsection |l fB|, i.e. 



n 



d=6,r=4 

[IV 



(?) 



d 3 k 



\k\<kp 



(2tt)3 2E k 



T, 



s=2,r=4 
fiv 



(q,k) 



(3.42) 



with the twist-4 spin-2 contribution to the forward scattering amplitude between a nucleon 
and the respective isospin current (|2.2|), 



ps =2,r=4/ 

- [IV 



(q,k) := i d 4 xe lllx (N(k)\Tj,(x)j u (0)\N(k)) 



s=2,t=4 



(3.43) 



With the general ansatz (cf. ( |3.20| )) 

T°f' T= \q, k) = ^ [q 4 S^(k) - q%q a S va (k) - q 2 [ q u q a S '^(k) + g^q 2 q a q p S a p{k)} 



+ B 2 {q^qv - 9[mvQ 2 ) q a q f3 S al3 (k) 
we get for the kinematical situation of DIS 

7jf =2 ' T=2DIS ( 9l k) = -q 2 (q ■ kfiB, d, u + B 2 e 
which has to be compared with eq. (2) from |^0| . We end up with 

B x = 1 (a 1 + -A 2 + -A 9 



[IV ) 



q 



and 



B- ? 



1 C- A 2 - l A ' 
q s U 8 



(3.44) 



(3.45) 



(3.46) 



(3.47) 



Based on a flavor decomposition new parameters KL are introduced in f4(| in terms of which 



the coefficients A J can be expressed. We refer to |^0[ for details and only give the final result 
for our isospin averaged coefficients A 3 : 



A 1 = -[K l u + K d - {l±l)K l ud }, 



(3.48) 



A 2 = -(K 2 U + K 2 ) 



(3.49) 



1 

3' 



A' = -IK' U + Kf 



d) i 



(3.50) 



where the + sign in ( |3.48| ) corresponds to the p and the — sign to the u meson. The 
parameters KL are given by the expectation values of the twist-4 condensates (|3.36| - |3.38" 
with respect to a single proton state p (cf. ( |3.36H3.38| , |3.41| )): 
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(p(k)\STg 2 u W X a u (u lvl5 X a u + d lv j 5 X a d)\p{k)) =: K x u S^{k) 
(j>(k)\STg 2 Ti ltl \ a u (u lu \ a u + d lv \ a d)\p{k)) =: K 2 u S^(k) , 



(3.51) 
(3.52) 



(p(k)\ST ■igu{D lt ,G va }'f'y 5 u\p(k)) =: K^S^k) (3.53) 

and respective definitions for K d , j = 1,2, g. Furthermore, 

(p(k)\STg 2 u W \ a u d lul ,\ a d\p(k)) =: K l ad S, u (k) . (3.54) 

Since the flavor structure of K d ' 2,g and K x ' 2,9 are governed by the d quark and the u quark, 
respectively, it seems reasonable to assume that the ratio is always the same 

(3.55) 



Ki Kl Kl - p - 



Within that assumption, the ratio (3 and the quantities B\ and B2 for p and uj can be 
determined from the DIS data. As expected from the valence quark decomposition of the 
proton one finds 



P « 0.5 



within a small error. Furthermore we get 

1 



B 1 = -[18(c l -c 2 )-(l±l)K 1 n 



ad I 



and 



B, 



;Q c 3 



(3.56) 



(3.57) 



(3.58) 



where the constants q, 2, 3, as well as K\ d can be found in table |. 

In total, the contribution of the twist-4 spin-2 operators is given by (|3.42j , l3.44l , l3.57t[338"|) . 
We note that a difference between the p and the 10 meson within the OPE up to dimension- 
ality d=6 only shows up for the twist-4 condensates and is expressed here in terms of the 
quantity K x ud . 

To summarize, we have presented in this section the operator product expansion of the 
current- current correlator ( |2.1D including condensates up to dimensionality d=6. The general 
form for the transverse and longitudinal part of the current-current correlator which enter 
( p. 15] ) is given by 



■in' 



1 + 



7T 



111 



'Q 2 



+ E 



-T,L 



Q 



2/i 



(3.59) 



where the coefficients have to be deduced from the various contributions discussed in 
this section. We note that contrary to the vacuum case PD|P^| and to the in-medium case 



with vanishing three-momentum q [25| , ^3 , ^2 , 26| , 28l , ^5 , ^0| for the general case q 7^ there 
are not only l/Q 4 and 1/Q 6 terms in ( |3.59| ) but also higher order terms, even when we 
restricting the OPE to d<6 condensates. Appropriate powers of if in the numerator serve to 
achieve the correct overall dimension. We will come back to that point in section |V). 
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IV. QCD SUM RULE 



In the last section we have calculated the l.h.s. of ( |2.15| ) within the operator product 
expansion and some additional assumptions. We postpone the discussion of these assump- 
tions to section [VT| and present here some general ideas about the calculation of the r.h.s. of 
( [2.15| ) and about the use of this equation. 

If one has a model at hand which yields the current-current correlator for arbitrary 
positive energy and arbitrary three-momentum, one could directly use (|2.15|) to judge the 
reliability of this model. In practice, however, the situation is such that one might have a 
model for the lowest hadronic resonance in the respective isospin channel, i.e. for p and lo, 
respectively, but one usually has no model which remains valid for arbitrary high energies. 
In the dispersion integral of (|2.15|) higher lying resonances are suppressed, but only by a 
factor 1/s 3 . Clearly, it is desirable to achieve a larger suppression of the part of the hadronic 
spectral function on which one has less access. To this aim, a Borel transformation |20|.|37 



can be applied to ( |2.15|) . For an arbitrary function f(Q ) the Borel transformation is defined 

as 



f(Q 2 ) 4 f(M 



with 



B := lim i 

Q 2 ^oo,N^oo T(N) 
Q 2 /iV = :M 2 =fixcd 



JV 



-Q 



2\N 



dQ 2 



where M is the so-called Borel mass. 

We will apply the Borel transformation to (cf. Q2.15 . 3.5S )) 

c n T . L (q 2 ) 




Q 



2n 



Q 2 



ct,l(<? 2 ) + — Ida 

71 J 



mill 



HAD 



T,L\s,q- 



[s + Q 2 -ie)(s + ie') 2 



Therefore we need to know the Borel transforms of f{Q 2 ) = {Q 2 + s) 13 and f{Q 2 ) 
From the definition fl4.2|) it is easy to derive [37| 

1 1 



f(Q) = (Q 2 + s) 



f(M 2 



r(/3) m 2/3 



-s/M 2 



and 



f(Q 2 )=lnQ 2 =► f(M 2 ) 
Applying the Borel transformation to ( |4.3|) we get 
1 ^ a s \ ^1 c^ L {q 2 ) 



-1 



8tt 2 



+ vr J + ^ Tin) 



n HAD^n ^ 



\0, q' 
M 2 



7lM 2 



M 2n 

ds 



ImU^is^q 2 ) 
s + ie 



s/M 2 



(4.1) 



(4.2) 



(4.3) 
InQ 2 . 

(4.4) 
(4.5) 



(4.6) 
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It is useful to write the r.h.s. of the last equation in a form, where it is more obvious that 
this expression is actually real valued. To this aim we split l/(s + it) into a principal value 
and a 5-function: 

' " m8(s). (4.7) 



s + ie s 2 + e 2 
Using this decomposition we find the QCD sum rule 



8vr 2 V * J n r W M 2 
Ren?AD( 0) g'2 ) 1 



J dsImn*^(s,q 2 )^— 2 e-*/ M2 . (4.8) 



M 2 ttM 2 J '' L v V + e 

We observe that higher resonance states are now exponentially suppressed. Additionally 
we find a 1/s suppression. The latter is due to the fact that we have applied the Borel 
transformation to l/Q 2 times equation (|2.15|) instead of directly applying it to (|2.15|) . On 



the one hand, such an additional suppression factor is desirable. On the other hand, we 
have to pay a price for this, namely that the subtraction constant n^ D (0, q 2 ) has not 
dropped out in contrast to the other subtraction constant ct,l(<? 2 )- Had we applied the 
Borel transformation to l/Q 4 times equation ( p. 15 ), the latter would also have survived. 



This is of course easy to understand from the point of view of subtracted dispersion relations: 
the suppression of high energy contributions has to be compensated for by a more detailed 
knowledge of the function at the subtraction point. We note that it is easy to get from 



also the direct Borel transformation of ( 2.15 ) without the l/Q factor. We simply have to 



multiply (^4.8| ) by (— M 2 ) and differentiate with respect to M 2 afterwards. Using this recipe 
the subtraction constant IT^ D (0, q 2 ) obviously would drop out. Also the 1/s suppression 
in the integral of (|4.8| ) would disappear. 

Having achieved a reasonable suppression of the energy region above the lowest lying 
resonance the integral in ( |4.8|) is no longer sensitive to the details of the hadronic spec- 
tral function in that region. For high energies the quark structure of the current-current 
correlator is resolved. QCD perturbation theory becomes applicable yielding 



ImIT™( S , q 2 ) = ( 1 + — ) for large s. (4.9) 

These considerations suggest the ansatz 

ImIT»l D ( S , q 2 ) = G(s - a) ImIT* E L s ( S , q 2 ) + ®{s - s ) Zf (l + 2i) , (4 .1 ) 

where So denotes the threshold between the low energy region described by a spectral func- 
tion for the lowest lying resonance, ImIl RES , and the high energy region described by a 
continuum calculated from perturbative QCD. Of course, the high energy behavior given in 
( OH ) is only an approximation on the true spectral function for the current- current cor- 



relator. Also the rapid cross-over in ( [4.10|) from the resonance to the continuum region is 



not realistic. However, exactly here the suppression factors discussed above should become 
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effective making a more detailed description of the cross-over and the high energy region 
insignificant. 

The price we have to pay for the simple decomposition ( |4.10| ) is the appearance of a new 
parameter Sq, the continuum threshold, which in general depends on the three-momentum 
q and on the nuclear density. We will elaborate later on the determination of so- 

Inserting (|4.10| ) into flOD yields 



8tt 2 



1 + — 

7T 



-s (q 2 )/M' 



E 



"T,L 



(? 2 ) 



Y{n) M 2n 



ReU^f{0,q' 
M 2 



s {q 2 ) 



nM 2 



dslmU^f(s } q 2 ) 



S 2 + 6< 



s/M 2 



(4.11) 



Obviously, the exponential suppression in ( [4.11| ) works only, if the Borel mass M is not 
too large. On the other hand, the OPE expression on the l.h.s. of ( |4.11|) gives a reliable 
prescription, if M is not too small, since we have neglected higher order condensates which 
are accompanied by higher orders in 1/M 2 . At best, the sum rule ( |4.11|) is valid inside of a 
Borel window 



M ■ < M < M , 

min — — max ' 



(4.12) 



where M^ in has to be determined such that the neglected condensates do not spoil the 
validity of the l.h.s. of ( |4.11| ), while M^ ax has to be determined such that the suppression of 



the details in the high energy structure of the current-current correlator becomes effective. 
Of course, it is not clear a priori, if such a Borel window exists at all. It might happen that 
M min is larger than M max . In this worst case, the sum rule ( |4.11| ) would be useless. 
The strategy to determine the Borel window is as follows f22|| : 
- For M min we require that for this Borel mass the absolute value of the contribution of the 
d=6 condensates is a certain percentage p of the total absolute value of the l.h.s. of ((4.8|). 
Since the d=6 condensates have the highest order in mass which is taken into account, one 
might expect that the relative contribution of the neglected condensates is much less than 



p. Following [22 1 we take p 



10%, i.e., 



n, d= 
Z T,L 



(M- 



min) 



= 0.1 








87T 2 V 7T J 



(q' 



r r(n) (M; 



min) 



(4.13) 



- For M max we require that for this Borel mass the absolute value of the continuum con- 
tribution to the integral in (|4.8| ) is a certain percentage p' of the total absolute value of the 
integral. Again we follow [B3] and take p' = 50%, i.e., 



J dse(s-s )ImIT^ D (s,g : 



s 2 + e 2 



0.5 



ds Indies, q' 



s 2 + e 2 



ax 



(4.14) 



By insertion of the decomposition ( [4.10| ) we get 
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J_ ( 1 + Si) g-off 2 )/^ = _ f dg i mU Rff Si f)-^— e ~ s / M ™ (4.15) 
8n \ n J J ' s 2 + e 2 

We note in passing that the sign of Imn^ D and therefore also the sign of Imllj; 1 ! 3 is always 
negative (cf. (|4.9| J2.12|) ). 



Obviously, the lower limit of the Borel window depends only on the condensates cal- 
culated in section |TTT|. In contrast to that, the upper limit depends on the choice for the 
continuum threshold So and on the hadronic model which yields E^f. In general, both 
limits may depend on the three-momentum q and on the nuclear density. 

Figure [l] shows the transverse component of the l.h.s. of (|4.8|) as a function of the Borel 
mass squared, M 2 , for various values of three- momentum squared, q 2 , and for p and uj 
meson. For the nuclear density we have chosen the nuclear saturation density of 0.17fm -3 . 
Figure ^ shows the same for the longitudinal component. On the left hand side both figures 
start with M 2 = M^ in as deduced from (|4.13|). Obviously, the difference between p and uj 



meson is only very small and vanishes with rising M 2 . The latter observation can be easily 
understood recalling that the only difference in the OPE's for p and uj comes from the twist-4 
spin-2 condensates which are suppressed at least by a factor 1/M 6 . Hence, the suppression 
becomes more effective with rising M 2 . Note that the small difference between p and uj 
does not necessarily mean that there is not much difference in their spectral functions. It 
only means that the integrated quantity given in ( |4.11| ) (to be rigorous: the r.h.s. of that 



equation) is nearly the same for both mesons. Nonetheless, the fact that the sum rule 
is nearly insensitive to the choice for the meson provides a strong constraint on a hadronic 
model which aims at a description of p and uj on the same footing, like e.g. ||15|| . 

We also observe that the dependence of the l.h.s. of the sum rule (|4.8|) on the three- 
momentum q is rather weak. This also constrains the hadronic model. Again, we stress that 
this does not mean that the dependence of the spectral function on the three-momentum is 
weak. 

Suppose now that one has a hadronic model at hand which yields at least the imaginary 
part Imllj^, (s, Q 2 ) f° r the respective isospin channel at finite nuclear density. Examples 
can be found in [|9^ |ll)| , |T7jj . Then, for given nuclear density and three-momentum q one can 
utilize the sum rule ( [4.11 ) and the results shown in figures [I] and ^ as a consistency check 



for the hadronic model in the following way (see also |30| 

Choose a continuum threshold so and a subtraction constant Rellj; 1 | s (0, q 2 ) 
Calculate the limits of the Borel window according to ( [4.13|) and (^4.15|) . 



Calculate the relative deviation r of the r.h.s. from the l.h.s. of the sum rule ( [4.11 ), 
averaged over the Borel window, i.e. schematically 

M 2 

max 

r= J d(M 2 ) |1 -r.h.s./l.h.s.|/AM 2 (4.16) 



with 
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(4.17) 



The l.h.s. function can be taken from figures |T| and ^. 

• Tune the "fit parameters" so and Relief (0, q 2 ) such that the deviation r becomes 
minimal. 

If this optimal r is reasonably small and the Borel window not too small, one might conclude 
that the considered hadronic model is in agreement with the QCD sum rule for the chosen 
nuclear density and three-momentum q. 

We close this section with some remarks on the respective size of the "fit parameters" sq 
and Rellj; 1 | s (0, q 2 ). Clearly, the more fit parameters we have the less restrictive is the sum 
rule for the hadronic model which should be checked. At least, it is therefore important to 
get an idea about the size and the possible influence of the fit parameters. 

In vacuum the continuum threshold So turns out to be about 1 — 1.6 GeV 2 pO|Jl5| , p2| . At 
least it has to be below the exited states of p and uj. Model calculations suggest that the 
threshold decreases with increasing density [|25|j28|,p~5|j3C 



Concerning the subtraction constant Relief (0, q 2 ) it is important to note that within 
the Fermi gas approximation it can be rigorously calculated for vanishing three-momentum 
| 26Ul5l| . Here it turned out that it is so small that it would not change the results drastically, 



if it is simply neglected. Unfortunately, the expectation that it could be neglected also 
for finite three- momentum is presumably not justified. If we use the full electromagnetic 



current in (|2.2 ) instead of a part of it with well-defined isospin, then within the Fermi gas 
approximation the transverse part of the subtraction constant would simply be the real 
part of the forward scattering amplitude T(q, k) of a (real) photon with momentum q on a 
nucleon with momentum k averaged over the Fermi sphere, i.e., 

Ren r( 0,,-)=4_/ ^^RertM^ReTKO), (4.18) 

\k\<kp 

where we have used the linear density approximation (cf. next section) to obtain the last 
expression. It is reasonable to assume that ReITj; ES (0, q 2 ) is of the same order of magnitude 
as Relly m (0, q 2 ). Thus, to get an idea about the size of the corresponding quantity in the 
sum rule ( [4.1 1| ) we have plotted in figure [| the quantity 

1 PN -ReT(g,0) (4.19) 



M 2 2m N 



as a function of the photon energy E~ — \q\ for nuclear saturation density and a typical value 
for the Borel mass, M = 1 GeV. We have used the model for photoproduction presented 
in |^3| to obtain the real part of the isospin averaged photon-nucleon forward scattering 
amplitude. Comparing the absolute sizes in figure |T| and |3| we find indeed that the subtraction 
constant might be negligible for q = 0, but not for arbitrary three-momentum. Especially, if 
we are interested in the dependence on the three-momentum we have to take the subtraction 
constant into account, since the variation of the curves in [I] with three-momentum is of the 
same order of magnitude as the quantity plotted in figure H. Concerning the longitudinal 
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part of the subtraction constant we cannot compare with photon-nucleon scattering, since 
there are no real longitudinal photons. Therefore, we refrain from presenting any estimates 
for this case. 

Of course, the most fortunate situation would be, if the considered hadronic model 
already provides a value for Rellf; 1 | s (0, q 2 ). Then, only the continuum threshold sq would 
remain as a fit parameter. 



V. LINEAR DENSITY APPROXIMATION 



Obviously, the contributions of the OPE presented in section |TT| are quite unillustrative, 
simply due to its complexity. To get more insight in the various contributions we restrict 
ourselves in this section to the parts which are at most linear in the nuclear density, i.e. cubic 
in the Fermi momentum. Recalling that we have evaluated all in-medium condensates using 
the Fermi gas approximation ( |3.7| ) we have neglected nucleon-nucleon correlations anyway 
which are quadratic in the density. Thus, the results presented in section [HI] are at best 



correct up to Fermi momentum to the power of five, i.e. up to o(p^ ). Therefore, we do not 
lose too much information, if we restrict ourselves here to the linear density approximation. 
To put it in physical terms: What we neglect further on is the Fermi motion of the nucleons. 
Anyway, for concrete calculations we can use the full results presented above. We note 
in passing that actually all the required integrals over the Fermi sphere can be calculated 
analytically. The results, however, are lengthy and unillustrative. Thus, for pedagogical 
reasons it is useful to discuss the linear density case, 



d 3 k 



\k\<kf 



(2tt) 3 2£, 



■(N(k)\0\N(k)) 



Pn 
2m 



N 



(N(0)\O\N(0)). 



(5.1) 



In this case, the coefficients c^ L (g 2 ) introduced in (|3.59|) are given by 



c t,l(<? 2 ) = c T d L 4 (V 2 ) + Ct\l b (<? 



n, d=6 i 



(5.2) 



with the contributions from the d=4 condensates to the transverse part, 



n=2,d=4 



2l(v G2 ) + m ^ )0 



+Pn 



(o) 

27 2 4 



A% +d (C? 



2,2 



3 \ 2 \G i s~iG /"tO \ 



(5.3) 



n=3, 



d=4 



{q 2 ) = -pnq 



2 m N 



AT 



(5.4) 



and to the longitudinal part, 



n=2,d=A 



(5.5) 
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c 



n=3, d=ir ->2\ 



PnQ' 



L.2 



2 ^L,2 



(5.6) 



and the corresponding contributions from the d=6 condensates, 



n=3, d- 



+Pn 



112 . 2 
= — ^Ka s 7r(qq) 

112 oat ,_ v 5m% r 
Ka s n(qq) - 



81 m, 
+ 



i 

9m N 
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A i v^2,a 2 X ' 4 ' 1 3 4 ^ 2 > 4 2 ,4y 



c 1 -c 2 - l -c s -^{l±l)Ki d 



(5.7) 



n=4,d=6/ ->2 



(r 



9m 



AT 



4 

9mjv 



AU+d(/~tq 
/i 4 ^2,4 



01 - ca - ics - i(l ± 



(5.8) 



-2pArg 4 m^ 



A u + d (ni _ r" 7 1 -I- -A G (C! G — r G \ 

/i 4 [Ly 2 A ^L,i) + l°2,4 °L,4J 



(5.9) 



and 



n=3, d= 



n=3, 



(5.10) 



n=4, d=6f ->2 



(q ) = pnq ' 



m 



N 



Au+d/fiq _ r-pq \ ZA u (n u — ^r G 



3mArC3 



(5.11) 



r 5 ' d=6 (? 2 )=2p^g 4 m^ 



9 

/iM+d^yg 1 ^ aGsiG 
^4 U L,4 "T 2^4 °L,4 



(5.12) 



All coefficients which are not given explicitly above vanish. As always, the ± sign which 
accompanies the K\ d term corresponds to p and u. 

We immediately observe that the scalar condensates contribute only to the three- 
momentum independent coefficients c?f~l' d ~ A and c^f' ~ 6 . These coefficients are identical 
for transverse and longitudinal part, since at q = we cannot distinguish between transverse 
and longitudinal directions |23[ . 

Obviously, there are contributions from numerous condensates to each of the coefficients 
presented above. On inspection of the parameters given in tables I and we can work out 
which condensates dominate which coefficient. Neglecting less important condensates we 
find 



n=2, d=i 



L 24 \ 7T /o 4 2 ' 



i+d 



(5.13) 



7i=3, d=4 
C'jp 



-PnQ' 



,m N 
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(5.14) 
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n=3,d=4s ->2\ 



n=3,d=6 n=3,d=6 

C T — C L 



, . I n=3,d=4/ ->2 



(? 2 )l 



112 2 112 a* 

— Ka s 7r{qq) - piv — Ka s ii{qq) 



(5.15) 
(5.16) 



ci-c 2 -ic 3 -^(l±l)^ d 



(5.17) 



c r 5 ' d=6 (g 2 ) « -2p N q A m\Al +d 



n=4, d=6 / ->2 



(5.18) 
(5.19) 



=5,d=6/- r 2> 



n=5, d=6 
T 



(<f 2 )l 



(5.20) 



Especially we have neglected all a s corrections to the twist-2 condensates, i.e. we have 
approximated Cf n by 1 and neglected all other C l rn (cf. ||35||). For vanishing three- momentum 
q the density dependent terms are dominated by the twist-2 spin-2 quark condensate and 
the scalar four-quark condensate. Concerning the q 2 -terms it is remarkable that the twist- 
4 condensates are equally important as the twist-2 spin-4 quark condensates. Of course, 
both are suppressed for large Borel masses as compared to the (d=4) terms. Thus, for 
the transverse part the twist-2 spin-2 quark condensate governs the g* 2 -terms. For the 
longitudinal direction the situation is more involved, since the (d=4) coefficient given in 
Q5.6D is quite small. Therefore, we have competing contributions from (|5.6j ) and (|5.19 ). 

We stress that in principle it is not necessary to perform the approximations which have 
led from (|5.3| - |5.12|) to (|5.13| - |5.20D . Of course, one can use the exact expressions for the 
coefficients. The purpose here was only to figure out the condensates which have the most 
influence on the coefficients. 

Since we have non- vanishing coefficients up to n=5 we find contributions up to o(l/Q 10 ). 
One may suspect that it is inconsistent to keep terms of order o(l/Q s ) and higher, since 
we have neglected (d—8) condensates which would contribute at o(l/Q s ). However, this 
is misleading, since the dependence on the three-momentum q is different in both cases. 
Schematically the neglected higher order condensates would contribute as 



(d = 8) condensate 



Q 



Q 2 



(5.21) 



where # denotes arbitrary dimensionless numbers which do not depend on q. Thus, e.g. the 
g 2 -terms of the neglected condensates are actually o(l/<5 10 ), while the corresponding terms 
of the condensates taken into account are o(l/Q s ). 

To work out the dependence of the current- current correlator on the three-momentum q 
more explicitly we study ( |4.8j ) in the vicinity of q = 0. For vanishing three-momentum we 
find 
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1 / 1 

7T 

rHAD 



^,n=2,d=4 ^,n=3,d=6 
+ TT^~ + 



M 4 



Ren HAU (0,0) 
M 2 



irM 2 



2M 6 
cMmIT HAD (s, 0) 



s 2 + e 2 



-s/M 2 



(5.22) 



We have skipped the label T, L, since there are no distinct directions at vanishing three- 
momentum. Next we differentiate ( |4.8|) with respect to q 2 and put q = afterwards. This 
yields 



in=3, d=4 in=A, d=6 
a T,L , a T,L 



2M 6 
d 



d(q- 



6M 8 



RenKi D (0,g 



V 



M 2 



irM 2 



rf S Imn«l D ( S ,g 2 )^^ 
s z + e z 



-s/M 2 



(5.23) 



g=0 



with 



d: 



n=3,d=4 U T,L 



^n=3,d=4/^2 



(q- 



q- 



In the same way we find by differentiating twice: 



d 



ii=5, d=6 



d 2 



24M 10 d(q 2 ) 2 



with 



V 



Ren^l u (0,g : 
M 2 



n=4, d= 



n=4,d=6 _ C T,L 



ttM 2 



dslmU^is, q- 



(5.24) 



s 2 + e 2 



-s/M 2 



g=0 



(5.25) 



m=5. 
a T,L 



=5,d=6 



(5.26) 



Let us discuss now the range of validity for the new sum rules ( |5.23j ) and ( |5.25| ) . As pointed 



out in section |T| it is crucial to find a non-vanishing Borel window where the validity of 
the sum rule is guaranteed inside of this window. To find the lower limit of this window we 
have compared the contribution of the highest order condensates to the l.h.s. of ( |4.8| ) with 
the total result (cf. ( ^.13| )). In ( |5.22| ) we have four different orders in 1/M 2 , namely zeroth 
to third order.0 Thus, it is no problem to compare the third order contribution to the total 
result. In ( |5.23|) we are left with third and fourth order in 1/M 2 , only. Thus, the number 
of orders we have access on is already diminished. This leads to a lower limit of the Borel 
3 GeV 2 for the transverse component of sum rule (|5.23|) which is already 



window of M^ in pa 



5 Note that the first order term vanishes. Strictly speaking it is proportional to the light current 
quark mass squared which is negligibly small. 
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much higher than the one for sum rule (|5.22 ): M^ in rs 0.6 GeV 2 . For the longitudinal 



component of sum rule ( |5.231 ) we even find M^ in m 10 GeV 2 . As already discussed above the 
q 2 -part of the (d=4) contribution to the longitudinal part ( |5.6| ) is quite small. Therefore, 
only for very large values of the Borel mass the (d=4) contribution can overwhelm the (d=Q) 
contribution. Both for the longitudinal and the transverse part we find that the respective 
lower limit of the Borel window for sum rule (|5.23| ) is much higher than the one for ( |5.22|) . 
If the upper limit of the Borel window does not rise in the same way, the sum rule ( |5.23| ) 
would not be as useful as (|5.22 ). To determine the upper limit of the Borel window we 
would need, of course, a hadronic model. For the most simple case, the approximation of 
the spectral function by a 5-function, it was found in |29[] that the upper limit of the Borel 
window also rises. Thus, also the sum rule ( |5.23| ) might be useful as a consistency check for 
hadronic models. We note, however, that the definition of the Borel window in [29|] differs 
from ours. For (|5.25 ) the situation is even worse. There we have only access on one order in 
1/M 2 . Thus, we cannot determine a lower limit for the Borel window. The sum rule ( p.25| ) 
is therefore useless. 

We stress again that the approximations performed in this section are not mandatory. 
The purpose of these approximations was to obtain more qualitative insight in the impor- 
tance of the various contributions and in the dependence on the three-momentum q. To 
check the consistency of a hadronic model with the QCD sum rule ( 4.11 ) the OPE coeffi- 
cients should be deduced from the formulae presented in section |Hl| . Only if the hadronic 
model is also restricted to the linear density case, a direct comparison with the simplified 
expressions would be appropriate. 



VI. DISCUSSION OF THE ASSUMPTIONS 



In this section we will discuss the various assumptions that have led to the condensate 



contributions calculated in section III 



The basic assumption we have made for the evaluation of in-medium condensates is the 



Fermi gas approximation (3/7 ). Clearly, this approximation is only valid for not too high 
densities. Of course, it would be of interest to quantify this statement. Unfortunately this 
is hard to do within the OPE approach, since it is not clear how to calculate expectation 
values with respect to multi-nucleon states for arbitrary nuclear densities. An idea about the 
importance of multi-nucleon states might be obtained from a comparison of the moments 
of parton distributions ( |3.26j , |3T2"TD as deduced from DIS experiments with nucleons on the 
one hand side and with nuclei on the other. This might be a direction of future studies. 
Concerning the parton distributions in nuclei we refer to MM and references therein. 

Within the framework of a hadronic model it is tedious but possible to approximately 
take into account interactions of the vector meson with more than one nucleon [p|-|T4],17l. In 



TTfl it was found that already at nuclear saturation density it is important to account for 



such processes. On the other hand, presumably in every hadronic model one can distinguish 
between single and multi-nucleon interactions. Therefore, it should be possible to compare 
the OPE calculation with the hadronic model treated in a "single nucleon mode". This 
comparison was e.g. performed in [T5[] for vanishing three-momentum q. If in the framework 
of a hadronic model it turns out that interactions with more than one nucleon are important 
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for the nuclear density under consideration, then the QCD sum rule cannot serve to check 
the consistency of the whole hadronic model but only of its restriction to scattering processes 
of the vector meson with one nucleon from the Fermi sphere. 

After these general considerations about the validity of the calculation of in-medium 
condensates we turn now to the discussion of the accuracy in the determination of the 
different types of condensates. 

Concerning the scalar condensates of subsection |111 A| we have neglected there terms 
which are quadratic in the light current quark masses m q as well as possible differences in 
the condensates of up and down quarks. In view of the fact that the light current quark 
masses are about 6 MeV, while the Borel window for the QCD sum rule ( |4.8| ) starts at about 
M^ in ps 0.6 GeV 2 (cf. figures [l|, ||), the neglect of m 2 q terms is very well justified. Also, a 
possible difference between up and down quark condensate is presumably small and anyway 
hard to disentangle in view of the uncertainties in the determination of the (average) light 
quark condensate (see e.g. discussions in j3U].f22|] ) . In addition, the contribution from the 
two-quark condensate to ( |3.6|) is much smaller than the one from the gluon condensate. The 
largest uncertainty lies in the evaluation of the four-quark condensate, i.e. in the value for 
k. First, even the vacuum value is still under discussion. Second, it might very well be 
that the value for k varies with nuclear density. Concerning the first problem, it is useful 
to choose a hadronic model which describes the data for e + e~ — > hadrons reasonably well 
and utilize the sum rule (|4.11|) for the vacuum case to determine k. This was performed in 



15] and we therefore have taken the condensate values given there. For the second problem 



we have no solution to offer so far. Without any better knowledge we use the vacuum value 
for k also at finite density. To get rid of this uncertainty one can differentiate the sum rule 
(4.8) with respect to the three-momentum squared q 2 . In this way, all scalar condensates 



drop out, since they do not yield <f 2 -dependent contributions to the OPE side of the sum 
rule. As already discussed in section [V] this differentiated sum rule results in a lower limit 
for the Borel window which is much higher than the one for the original sum rule. This 
might diminish or even close the Borel window, i.e. the thus obtained sum rule might be 
less reliable or even useless. This clearly depends on the explicit hadronic model under 
consideration. 

The twist-2 condensates discussed in subsections [MB , |III Q are the best known con- 
tributions, as soon as one accepts the Fermi gas approximation discussed above. We even 
can use the results of DIS to get an idea about the neglected higher order condensates (see 
below) . 

The twist-4 spin-2 condensates can in principle also be deduced from DIS data. The 
uncertainties in their extraction are, however, quite large. We have adopted the analysis 
of [|0|. There, condensates depending on the light current quark masses, strange quark 
contributions and dependences on the renormalization scale were neglected. These errors are 
presumably smaller than the uncertainties in the extraction of these condensates from DIS 
data. In general, the contributions of the twist-4 spin-2 condensates are small as compared 
to the twist-2 spin-2 contributions. For the <f 2 -contributions to the longitudinal direction, 
however, the twist-2 spin-2 contribution is suppressed by a s . There, the twist-4 spin-2 



condensates cannot be disregarded [^9[ . 

Of course, one is forced to somewhere cut off the OPE used here as an expansion in the 
dimensionality of the condensates. We have neglected all condensates of dimensionality 8 
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or higher. To get an idea about the size of the neglected condensates we calculate now the 
twist-2 spin-6 contribution to the current-current correlator, since we have access on that 
quantity utilizing the DIS results. The calculation proceeds along the lines described in 
subsections |III B| , [III Q We use the Fermi gas approximation 

(eW«*>«4 / ^^ k (N(k)\0, VKXix \N(k)) (6.1) 

\k\<k P 

and the decomposition 
(N(k)\0^ x \N(k)) 

~ k^kuk^kxk^ - —{k^k v kJz\g^ x m 2 N + 14 permutations) 

+ j^(k^k v g K x9^x m N + 44 permutations) - — (g^g K xg^ x m N + 14 permutations) 
=: S^ x (k) . (6.2) 

The relevant operators are 

ST t(u llM D u D K D x D ( D x u + d llM D u D K D x D^D x d) (6.3) 

and 

STG a ^D v D K D x D i G ax . (6.4) 

From the DIS calculations we can deduce the following contribution to the current-current 
correlator ( |2.1| ) which was neglected in ( |3.3| ): 

ifci<fc F 

with the forward scattering amplitude 

64 

T^ s =^= 2 (q, k) = - — (q^q^ - g«q 2 q u qP - g u a q 2 q^ + gfg/q') ffV W 

x (cl,Al +d + |c£^) 

--^(v^v ~ 9^q 2 ) q a q p q 1 q 5 q e q^S a p~ fS< (k) 

x (Vie - Cl,)Al +d + |(Cg 8 - C^) • (6-6) 

With this at hand we can calculate the ratio between the twist-2 spin-6 contribution to the 
l.h.s. of the sum rule 04.81) and the total value for this l.h.s. as calculated in section ITTTl We 



have calculated that ratio in figures f|, |5| for transverse and longitudinal directions and for 
p and us meson. Obviously, the obtained ratios are very small justifying at least the neglect 
of twist-2 spin-6 condensates and suggesting that all higher dimensional condensates are 
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reasonably suppressed. Of course, for a thorough discussion of the error made by neglecting 
higher dimensional condensates we also would need to know the other condensates beside 
the twist-2 spin-6 condensates. Since we do not know e.g. the scalar (d—8) condensates 
etc. the error analysis presented here is only a first guess. 

As long as we do not know the actual values of the neglected condensates (or at least 
an upper limit for them) we cannot present a rigorous proof that the OPE approach works, 
i.e. that the truncated series yields a reliable value for the current-current correlator in the 
region of interest. Indeed, in it was doubted that the QCD sum rule approach provides 



any reliable information about the medium modifications of vector mesons (see also [f4"5|,[4"6]]). 
Two arguments were given there to support these doubts: The first, qualitative argument 
concerns the connection between the mass shift of a vector meson and the forward scattering 
amplitude of this vector meson with a nucleon. It was argued there that the forward scat- 
tering amplitude and hence also the mass shift is a long distance property, while the OPE is 
only capable of short distance properties. We think that this argument is misleading, since 
the OPE always concerns the interplay of long and short distance properties, as already 
pointed out in section pij . Actually, in the same oversimplifying way one might argue just 
the other way round: The in-medium p mass is still large, i.e. a short distance property, and 
therefore can be described by the OPE approach. This shows that one needs more quan- 
titative arguments to check the validity of the sum rule approach. A useful selfconsistency 
check within the Borel sum rule method is the Borel stability analysis described in section 
[V : A breakdown of the sum rule might be observed in a small or even vanishing Borel win- 
dow. Indeed, this stability analysis was the key point to resolve the question, whether mass 
shift and/or forward scattering amplitude can be extracted within the traditional sum rule 
approach utilizing the narrow width approximation. We refer to 1 47, 26, 28, 48] for details. In 
general, the Borel window can only be determined after specifying a hadronic model. There- 
fore, we do not discuss this point here any further. The preceding discussion has clearly 
emphasized the necessity to perform a Borel stability analysis. 

A second, formal argument has been raised in against the applicability of the OPE 
approach to vector mesons in nuclear medium: It was claimed there that the OPE turns out 
to be an expansion in the nucleon mass m N over the invariant mass \J~Q 1 . After Borelization 
this would turn into an expansion in mjv over the Borel mass M. If the latter is assumed to 
be of the order of the p meson mass, one would get an expansion parameter m^jm p which is 
obviously not small. Therefore, it was argued in |1| that the truncation of the OPE at the 
d=6 condensates is not appropriate. Indeed, concerning the twist-2 condensates the state- 
ment is true that the used OPE is an expansion in /M. This can most easily be discussed 
within the linear density approximation of section |V]. E.g. for vanishing three-momentum 
the twist-2 spin-2 condensates contribute with a term proportional to pN m N /M 4 (cf. flOp) 
and the twist-2 spin-4 condensates with a term proportional to p^m^/M 6 (cf. (|5.7j )). In 
general, the twist-2 spin-s condensates yield a contribution ~ p^m^ 1 /M s+2 . Since the 
nucleon mass is large this expansion might break down for the Borel masses of interest 
(typically of the order of 1 GeV). However, one should not discuss the convergence of a series 
without looking at its coefficients. The twist-2 spin-s contribution is accompanied by the 
s-th moment of the parton distributions. Inspecting the last two lines of table [TI| we find 
that these moments become very small with increasing s. Indeed, we have already discussed 
above that the twist-2 spin-6 condensates do not spoil the truncation of the OPE — in 
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spite of the fact that they are proportional to pN?n%/M 8 . This shows that also the second 



argument raised in [16| against the QCD sum rule approach is oversimplified. 



VII. SUMMARY 



This work was motivated by the finding that the QCD sum rule approach provides 
no model independent prediction about a possible mass shift of vector mesons in nuclear 
medium [jl5| , pD|| . In |]30| we have discussed at length that the sum rule restricts the p meson 
only to a rather wide area in the (mass, width) plane without making any further statements 
about the specific properties of the p meson. Only within additional model assumptions the 
behavior of the p meson in nuclear matter can be further specified. E.g., if one assumes 
that the width of the p meson is not increased, then the sum rule predicts a p mass which 
decreases with increasing nuclear density. However, it is also possible to assume instead 
that the p mass is not shifted. In this case the sum rule suggests an increasing width of the 
spectral function of the p meson. 

This, however, does not mean that the sum rule approach is useless: We have presented 
here a QCD sum rule for p and u mesons propagating with arbitrary three-momentum 
through nuclear matter at vanishing temperature. This sum rule provides a non-trivial 
consistency check for hadronic models which describe that propagation. At least as long as 
different hadronic models cannot be judged unambiguously by experiments such consistency 
checks are important to confirm or rule out hadronic models. 

The main formula is given in (|4.11 ). The OPE coefficients which appear on the 
l.h.s. of this formula are defined via n°£ E in ( |3.59|) . In view of their complexity we have not 
given the explicit formulae for L . However, they can be easily deduced in the following 
way from the equations presented in section HI]: Transverse and longitudinal part of 
are obtained from the respective last expression of (|2.7|) and (|2.8|). The current-current 



correlator with the full Lorentz structure, is decomposed in (|3. 



The scalar contribution is given in (3^) where _R scalar can be read off from (3J3). The 
expectation values showing up there are decomposed in vacuum and medium dependent 
expectation values in (|3.9| ) using the Fermi gas approximation. The vacuum expectation 
values are listed in table |. Finally, the medium dependent parts of the scalar condensates 
are connected in ( |3.10| , |3.1I| ) with parameters also listed in table |. 



The contribution from the twist-2 spin-2 condensates is given in (|3~28| , F29| ) with the 
traceless tensor S a p defined in ( |3.13| ). In the same way the contribution from the twist-2 
spin-4 condensates is given in (|3.34| , |3.35D with the traceless tensor S a p^& defined in Q3.31|) . 
The required values for the moments of the parton distributions A l n and the coefficients C l r n 
are collected in table 0. 

The contribution from the twist-4 spin-2 condensates is given in ( |3~4^ , CT ) where the 
coefficients -81,2 are connected in ( |3.57| , |3.58D to quantities listed in table |. 

In this way, the OPE coefficients can be easily calculated. As discussed in the last section 
the QCD sum rule ( |4.11| ) can be used to check the consistency of a hadronic model, provided 
that in the latter the medium is also described by the Fermi gas approximation. Going one 
step further by neglecting the Fermi motion of the nucleons both in the hadronic model 
under consideration and in the calculation of the OPE coefficients one might also utilize the 
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sum rule in this linear density approximation. For this case the OPE coefficients c^ L are 
explicitly given in (|5.2| - 5.12| ). 

By inspecting the QCD sum rule ( [4.11| ) we observe that the hadronic model which should 
be checked has to yield the current- current correlator for invariant masses in the region (— <? 2 ) 
to sq. The lower limit refers to vanishing energy. For non-vanishing three-momentum q this 
means that we need information not only about the time like, but also about the space 
like region. For small three-momenta the space like region is dominated by the coupling 



of the respective vector meson to nucleon-hole states p6[ . For higher three-momenta also 
resonance-hole loops come into play in the space like region (see figure 3 in |HJ). Thus, at 
finite nuclear density there are important structures in the spectral function of the vector 
mesons also in the space like region. This is in contrast to the vacuum case where there is 
no structure below the two-pion (three-pion) threshold for the p {uj) meson. 

Qualitatively, we have found that our sum rules are not very sensitive to the difference 
between p and uj meson as well as to a variation in the three-momentum of the vector meson 
with respect to the nuclear medium. This, however, does not a priori mean that the current- 
current correlator for different isospin channels and for different three-momenta should be 
more or less the same. In the sum rule only an integral over this correlator enters which 
might be the same e.g. for p and uj mesons, even if the respective correlators themselves are 
different. Thus, on this qualitative level the sum rule approach does not rule out hadronic 
models which predict a different behavior of vector mesons with different three-momenta, 



like e.g. ||13| , p^ , p!6H Tq] . A quantitative analysis of these models is beyond the scope of this 
paper. 

We believe that the QCD sum rule presented here provides an interesting and non-trivial 
consistency check for hadronic models which describe vector mesons in nuclear matter. We 
have tried to present the derivation of the sum rule in great detail to make it possible also 
for non-experts in OPE to utilize the sum rule for a consistency check of their hadronic 
models. 
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TABLES 





0.36 


a N 


[GeV] 


0.045 


Cl 


[GeV 2 ] 


0.005 


m q [GeV] 


0.006 


(0) 

my 


[GeV] 


0.75 


C2 


[GeV 2 ] 
[GeV 2 ] 
[GeV 2 ] 


0.011 


(qq) [GeV 3 ] 


-0.0156 










0.035 


((a s /7T)G 2 ) [GeV 4 ] 


0.012 










-0.088 


K 


2.36 















TABLE I. Parameters used in the calculation of the OPE contributions to the current-current 
correlator. See main text for details. 



^2,2 

^2,2 
Cl.2 



1.013 
0.050 
-0.042 
0.057 



°2,4 
^2,4 



1.171 
0.030 
-0.063 
0.023 



°2,6 
^2,6 



1.316 
0.022 
-0.060 
0.012 



u+d 



AU+d 

^6 



A, 



1.12 
0.83 



4 M+ 



0.11 
0.04 



0.03 
0.01 



TABLE II. Relevant coefficient functions Cf taken from [35| and moments of parton distri- 
butions A J n calculated from p| for /i 2 = 1 GeV 2 . 
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FIGURES 
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FIG. 1. Transverse part of l.h.s. of ( |4.§| ) as a function of the Borel mass squared, M 2 , for 
three-momenta \q\ = 0, 0.5, 1 GeV (top to bottom) and for p (full lines) and uj mesons (dashed 
lines). 
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FIG. 2. Longitudinal part of l.h.s. of ( |4.E| ) as a function of the Borel mass squared, M 2 , for 
three-momenta \q\ = 0, 0.5, 1 GeV (bottom to top) and for p (full lines) and ui mesons (dashed 
lines). 
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FIG. 3. Real part of isospin averaged fN forward scattering amplitude as a function of the 
photon energy in the rest frame of the nucleon (rescaled with / '(2mArM 2 ), see main text for 
details). 
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FIG. 4. Relative error made by the neglect of twist-2 spin-6 condensates in the calculation of the 
transverse part of the l.h.s. of ( ^4.8[) as a function of the Borel mass squared, M 2 , for three-momenta 
\q\ = (upper line), 0.5 (lower line), 1 GeV (middle line) and for p (full lines) and co mesons (dashed 
lines). The only noticeable difference between p and oj meson appears in the slightly different lower 
limits of the Borel window, i.e. in the starting points of the curves on the left hand side. 
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FIG. 5. Same as figure || for longitudinal part. The curves refer to |g| = 0, 0.5, 1 GeV (top to 
bottom). 
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